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ABSTRACT 

With the commissioning of XMM-Newton came the discovery of 2.5 — 10 keV spectral com- 
plexity in some narrow-line Seyfert 1 galaxies (NLSl). This high-energy complexity can be 
manifested as sharp, spectral drops or gradual curvature in the spectrum. Models which are 
normally considered are ionised reflection and partial-covering. In this work, we define two 
samples of NLSl: a complex sample whose members exhibit high-energy complexity (C- 
sample), and a general sample of NLSl whose 2.5 — 10 keV spectra do not strongly deviate 
from a simple power law (S-sample). We than compare multiwavelength parameters of these 
two samples to determine if there are any distinguishing characteristics in the complex NLS 1 . 
Considering historical light curves of each object we find that the C-sample is representative 
of NLSl in a low X-ray flux state, whereas the members of the S-sample appear to be in 
a typical flux state. Moreover, from measurements of aox with contemporaneous UV/X-ray 
data, we find that the C-sample of NLS 1 appear X-ray weaker at the time of the observation. 
For two NLS 1 in the C-sample multi-epoch measurements of aox are available and suggest 
that aox approaches more normal values as the complexity between 2.5 — 10 keV diminishes. 
This implies that a source could transit from one sample to the other as its X-ray flux varies. 
Secondly, there are indications that the C-sample sources, on average, exhibit stronger optical 
Fe II emission, with the three most extreme (Fe ll/H/3> 1.8) Fe II emitters all displaying com- 
plexity in the 2.5 — lOkeV band. It is an intriguing possibility that we may be able to identify 
X-ray complex NLSl based on the extreme strength of the more easily observable optical 
Fe II emission. However, it is not clear if the possible connection between Fe II strength and 
spectral complexity is due to the Fe II producing mechanism or because strong Fe II emitters 
may exhibit the greatest variability and consequently more likely to be caught in an extreme 
(low) flux state. Based on the current analysis it we can not straightforwardly dismiss absorp- 
tion or reflection as the cause of the X-ray complexity; by considering the multiple UV/X-ray 
observations of IH 0707^95 (a C-sample member), we discuss a possible method of distin- 
guishing the two models provided further UV/X-ray observations. 
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1 INTRODUCTION 

The importance of narrow-line Seyfert 1 galaxies (NLSl) be- 
came apparent on the discovery of the 'primary eigenvector' (PCI; 
Eigenvector 1) in the Boroson & Green (1992) principle compo- 
nent analysis of PG quasars. PCI showed a strong anticorrelation 
between the strengths of [O III] and Fe II in the optical spectra, with 
NLSl as a class showing the strongest Fe II emission and weak- 
est [O III]. While the physical driver of PCI is still debated, there 
are strong indications that the fraction of the Eddington luminosity 
(L/LEdd) at which the object is emitting at is responsible, imply- 
ing that NLS 1 are relatively high accretion rate systems. Since there 
are no significant differences between NLS 1 and broad-line Seyfert 



1 galaxies (BLSl) in their X-ray, optical, or bolometric luminosi- 
ties (Grupe et al. 2004), it naturally follows from the condition of 
higher accretion rates that NLSl possess a lower mass black hole 
compared to BLS 1 of similar luminosity. Black hole mass measure- 
ments confirm lower mass black holes in NLSl (Wandel, Peterson 
& Malkan 1999; Peterson et al. 2000; Grupe & Mathur 2004). 

The strong Fe II emission in NLS 1 is prevalent in the UV (e.g. 
Laor et al. 1997; Constantin & Shields 2003), optical (e.g. Boroson 
& Green 1992), and infrared (e.g. Rodriguez- Ardila et al. 2002). 
The observed line widths and absence of forbidden emission sug- 
gests that Fe II is formed in the dense BLR, but photoionisation 
models cannot account for all of the Fe II emission. The 'Fe II dis- 
crepancy' remains unsolved, though models which consider non- 
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radiative lieating (probably due to shocks produced in outflows), 
with an overabundance of iron are promising (see Collin & Joly 
2000 for a review). 

Perhaps the most interesting characteristics of NLSl are man- 
ifested in the X-ray regime. The origin of the strong soft excess 
emission below about IkeV, known since Einstein observations 
(e.g. Puchnarewicz et al. 1992; see also Boiler, Brandt & Fink 
1996), is still debated (e.g. Mineshige et al. 2000; Gierlinski & 
Done 2004; Crummy et al. 2005); as is the nature of the ex- 
treme, rapid variability. With ASCA, came the discovery that the 
2 — 10 keV spectrum in NLSl also appeared steeper than in BLSl 
(Brandt, Mathur Sl Elvis 1997). The result could be understood 
as arising from significant Compton cooling of the accretion disc 
corona due to the strong soft X-ray excess found in most NLSl. 
The lower energy gain per scattering (smaller Compton y parame- 
ter) would yield steeper, hard X-ray slopes. 

With the high-sensitivity of XMM-Newton came what may 
be the most interesting discovery with regards to NLSl, and that 
was the presence of a sharp, spectral drop at i5 = 7keV in 
IH 0707-495 (Boiler et al. 2002). Since then, a number of drops 
at i? > 7keV (or more generally 'high-energy curvature') have 
been observed in several NLSl (e.g. Pounds et al. 2003, 2004; 
Longinotti et al. 2003; Boiler et al. 2003; Uttley et al. 2004; Reeves 
et al. 2004). The exact nature of this behaviour, which seems to 
be a characteristic of NLSl, is uncertain. Two models which ap- 
pear probable are partial-covering (e.g. Tanaka et al. 2004 and ref- 
erences within) and reflection (e.g. Fabian et al. 2002). 

In terms of partial-covering, the ~ 7 keV drop is produced by 
absorption of the continuum emission by a dense material, which 
only partly obscures the primary emitter This can, in principle, ex- 
plain the absence of other absorption features (e.g. intrinsic cold 
absorption, Fe L absorption, fluorescence emission), and if the ab- 
sorber is allowed to be in radial motion, it can possibly account for 
the various edge energies which are seen (e.g. Gallo et al. 2004a). 
It is important to realise that partial-covering does not describe the 
nature of the primary continuum source. Although a simple black- 
body plus power law continuum is often assumed, other physical 
descriptions of the primary continuum are not dismissed, nor do 
they discriminate against partial-covering. 

Reflection of the power law continuum source off the cold ac- 
cretion disc can also adequately describe the X-ray spectra of NLS 1 
(e.g. Fabian et al. 2004). In this case, the sharp drop at high ener- 
gies is the blue wing of a relativistically broadened iron line. In 
combination with light bending effects close to the black hole (e.g. 
Miniutti & Fabian 2004), the reflection model nicely describes the 
shape of the X-ray continuum and the principle of 'reflection dom- 
inated' spectra. 

It stands to reason that regardless of the correct model, the 
process may be ubiquitous in NLS 1 and probably in active galac- 
tic nuclei (AGN). By varying the prevalence of the physical pro- 
cess (e.g. the degree of absorption or amount of reflection) one can 
potentially describe the different types of X-ray spectra that are ob- 
served. Objects, such as IRAS 13224-3809 or IH 0707-495, mani- 
fest 'the process' significantly; thus exhibiting sharp and deep spec- 
tral drops. Other NLSl, like NGC 4051, demonstrate the process 
only moderately, only displaying gentle curvature over the hard X- 
ray band. In most other NLSl, the process is minimal and likely 
not detected. 

In this work we examine what, if anything, is unique about the 
NLSl which appear to possess high-energy complexity. Leighly & 
Moore (2004) began to address this issue when they examined how 
the UV properties of IRAS 13224-3809 and IH 0707-495 were re- 
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Figure 1. The null hypothesis of a power law plus Gaussian profile model 
to the intrinsic 2.5 — lOkeV spectrum of each NLSl in the sample plotted 
against 5 — 10 koV SN. The filled (black) circles mark the NLS 1 identified 
as complex (C-sample; see text for details). The open (red) squares are the 
general NLS 1 sample (S-sample). MCG-6-30-15 is shown as a blue star for 
illustrative purposes. 

markable compared to other NLSl. They found that IRAS 13224- 
3809 and IH 0707-495 possessed weaker emission lines, signifi- 
cant asymmetry in the C IV profile and steeper spectra. They con- 
cluded that many of the differences in the two NLS 1 could be ex- 
plained in terms of radiative line driven wind models (e.g. Proga, 
Stone & Kallman 2000; Leighly 2004). Here, we continue to probe 
this issue by investigating multiwavelength properties of a sample 
of NLS 1 , which includes a larger number of 'extreme' objects. 

The remainder of the paper is organised as follows. In the 
next section we define a complex sample of NLS 1 , which appear 
to exhibit high-energy complexity in their XMM-Newton (Jansen 
et al. 2001) spectra; and a general sample of NLSl, which do not 
show significant evidence of the complexity described here. In ad- 
dition, the X-ray data processing is described. The dependence of 
X-ray flux state in defining the samples is investigated in Section 3. 
In Section 4, the calculation of the optical-to-X-ray spectral in- 
dex from contemporaneous UV and X-ray data are described and 
presented. In Section 5 the multiwavelength parameters obtained 
from the literature (or estimated by us) are compared and our find- 
ings are highlighted. We discuss our results in Section|6|and give 
our summary in Section A value for the Hubble constant of 
Ho=70 km s^^ Mpc~^ and a standard flat cosmology with 
= 0.3 and Qa = 0.7 were adopted throughout. 



2 X-RAY DATA PROCESSING AND SAMPLE 
DEFINITION 

In defining a sample we considered all known NLS I that have been 
observed on-axis and in imaging mode with the EPIC pn camera 
(Striider et al. 2001) on XMM-Newton. The pn was selected to make 
use of what is most likely the highest signal-to-noise (SN), 2.5 — 
lOkeV rest-frame spectrum of these objects, currently available. 
At the time we commenced this study, there were more than 50 
observations of 37 NLSl -type objects (including quasars and low- 
luminosity sources). 

For each observation, the Observation Data Files (ODFs) were 
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Table 1. The X-ray data. The complete sample is divided into two sub-samples as described in the text. The information given is the: object name (col- 
umn 1); observation data (column 2) and XMM-Newton revolution (column 3); total good exposure for the pn (column 4); observed 2.5 — lOkeV flux 
(xlO~^^ ergcm"^ ) coiTected for Galactic absorption (column 5); intrinsic 2keV flux density (xlO"'^" ergcm"^ Hz~^ ) (column 6); null hy- 
pothesis from fitting the base model to the 2.5 — lOkeV data (column 7); 



General NLSl Sample (S-sample) 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


Source 


Date 


XMM-Newton pn exposure 2.5 


- 10 keV 


2keV 


NuU 




year.mm.dd 


rev. 


(ks) 


flux 


flux density 


Hypothesis 


Mrlc 7(^^^l 
iVllJv / UU 


7000 0^ 70 


0082 


26.0 


1 33 


19 08 


757 


iVliJv 


7000 07 OQ 


0107 


6.9 


51 


6.69 


418 


iVll Jv i L/'t't 


7007 07 73 


0480 


5.5 


57 


9.23 


578 




7000 C\h 1 7 


0096 


7.4 


1 86 


46.39 


992 


Mrk 896 


2001.11.15 


0355 


7.2 


0.31 


5.08 


0.339 


Mrk 335 


2000.12.25 


0192 


28 5 


1.20 


22.53 


0.731 


NGC 71 


200 1 1 1 27 


0361 


4 5 


0.02 


0.19 


0.890 


Mrk 493 


2003.01.16 


UJUo 


I J.J 


0.31 


5.95 


0.629 


IZw 1 


2002.06.22 


0464 


18 6 


0.68 


14.12 


0.182 


TonSlSQi 


2000.12.14 


0186 


20 6 


0.39 


9.00 


0.292 


PG 1448-H273 


2003.02.08 




1 O.J 


0.16 


3.71 


0.697 


UGC 11763 


2003.05.16 


0629 


26.2 


0.31 


3.46 


0.102 


RX 10323.2^931 


2003.08.16 


0675 


25.1 


0.11 


1.90 


0.103 


Mrk 478 1 


2003.01.07 


0564 


18.2 


0.15 


2.98 


0.972 


11 Zw 177^ 


2001.06.07 


0274 


8.7 


0.08 


2.09 


0.604 


IRAS 13349-1-2438 


2000.06.20 


0097 


30.4 


0.19 


3.08 


0.162 


PKS 0558-504^ 


2000.10.10 


0153 


8.0 


1.03 


20.75 


0.519 


PG 1115-H4072 


2002.05.17 


0446 


15.0 


0.10 


2.35 


0.746 


Nab 0205-H0242 


2002.07.23 


0480 


21.9 


0.13 


3.18 


0.320 


PG 2233-H1342 


2003.05.28 


0635 


7.5 


0.05 


1.15 


0.713 


E 1346-1-2662 


2003.01.13 


0567 


49.2 


0.005 


0.47 


0.464 






Complex 


NLSl Sample (C-sample) 








NGC 40511 


2002.11.22 


0541 


44.7 


0.59 


4.98 


< 0.01 


PG 1535-1-5472 


2002.11.03 


0531 


20.1 


0.12 


0.12 


< 0.01 


IH 0707^95 1 


2000.10.21 


0159 


38.1 


0.04 


0.39 


< 0.01 


IRAS 13224^3809 


2002.01.19 


0387 


54.3 


0.04 


0.84 


< 0.01 


PG 1211-H143 


2001.06.15 


0278 


49.5 


0.27 


3.24 


< 0.01 


PG 1402-H261 


2002.01.27 


0391 


9.1 


0.15 


3.30 


0.09 


PHL 1092^ 


2003.01.18 


0570 


19.8 


0.02 


0.31 


0.20 



FOOTNOTES: (1) Multiple XMM-Newton observations exist. (2) Adding an unresolved Gaussian profile was not an improvement over a simple power law fit. 



processed to produce calibrated event lists using the XMM-Newton 
Science Analysis System (XMM-SAS v6 . 1 . 0) and the most re- 
cent calibration files. Unwanted hot, dead, or flickering pixels were 
removed as were events due to electronic noise. Event energies 
were corrected for charge-transfer losses, and time-dependent EPIC 
response matrices were generated using the SAS tasks ARFGEN 
and RMFGEN. Light curves were extracted from these event lists to 
search for periods of high background flaring, which were then re- 
moved. The source plus background photons were extracted from 
a source-centred circular region with a radius of 35". The back- 
ground was selected from an off-source region with a larger radius 
and appropriately scaled to the source region. Single and double 
events were selected for the pn spectra. Pile-up was examined in all 
spectra and when corrections were necessary the inner pixels of the 
source extraction region were ignored. 

On processing of the spectral data, the sample was further re- 
fined by excluding observations in which: (i) the SN in the intrinsic 
5 — 10 keV range was < 1.2, or (ii) the intrinsic 5 — 9 keV spectrum 
was not source dominated. Finally, if more than one observation of 
the same object was available, the highest SN data were used. This 



filtering process resulted in the sample of 28 NLSl, which possess 
high-quality, intrinsic 2 — lOkeV spectra (Table0. 

The purpose of this study is to focus on the nature of NLS 1 
which exhibit complexity in their intrinsic ~ 2.5 — 10 keV spectra. 
To identify these objects, each spectra above 2.5 keV (rest-frame) 
was fitted with a baseline model composing of a power law mod- 
ified by Galactic absorption (Dickey & Lockman 1990) plus, if 
necessary, an unresolved (a < 150 eV) Gaussian profile between 
6.4 — 7 keV. The intention of this baseline model was not to deter- 
mine the 'best-fit' to the spectrum, but to identify the objects which, 
to first order, could not be explained by the most basic AGN model. 

Sources for which the null hypothesis of this baseline model 
was < 0.10 were marked as NLSl with high-energy complexity. 
We purposely adopted a low value for the null hypothesis to iden- 
tify the most extreme cases which, in turn, should make it easier 
to identify other extreme characteristic of the sample. Six NLSl 
(PG 1211-1-143, IRAS 13224-3809, IH 0707-495, PG 1535-^547, 
PG 1402-1-261 and NGC 4051) were categorised as having high- 
energy complexity. In addition, PHL 1092 was also included in this 
sub-sample because it exhibited a very flat, and likely unphysical, 
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2.5 — lOkeV power law (F ~ 1.55), even though the baseline 
model null hypothesis was 0.20. In total, seven objects make up 
the 'complex NLSl sample' (hereafter C-sample). The remaining 
twenty-one sources compose the 'general (simple) NLSl sample' 
(hereafter S-sample) (see Table0. 

Notably absent from the C-sample is the borderline NLSl, 
IRAS 13349+2438 (FWHM(H/3)^ 2800 kms"^ , but strong Fe II 
emission; Grupe et al. 2004), which clearly shows high-energy 
complexity (Longinotti et al. 2003), but does not satisfy the rather 
strict criteria established here. 

In FigureQthe null hypothesis is plotted against the 5—10 keV 
SN for the objects in the sample. As will be adopted throughout, the 
filled circles will mark the members of the C-sample and the open 
squares will identify the S-sample population. The distribution of 
SN is comparable for the C- and S-sample. Therefore, even though 
identifying X-ray spectral complexity obviously depends on SN, it 
does not appear that it will be a significant bias in selecting our sam- 
ple. For illustrative purposes the pn data for MCG-6-30-15 during 
revolution 302 are also included (blue star). We chose to include 
MCG-6-30-15 as it is perhaps the most obvious display of high- 
energy spectral complexity in a type 1 AGN. 



3 X-RAY FLUX DEPENDENCE ON SAMPLE 
DEFINITION 

We considered what effect X-ray variability could have on the sam- 
ple definition by examining the historic 2 — 10 keV light curve 
for each object. Specifically, we compared the XMM-Newton X- 
ray fluxes measured for each object with past measurements from 
other missions. 

ASCA fluxes were available either from Vaughan et al. (1999) 
or the TARTARUS archive for ten of the objects in the S-sample. 
Five objects also had BeppoSAX measurements (Comastri 2000) 
resulting in an X-ray flux measurement for eleven NLSl in the 
S-sample (four having BeppoSAX and ASCA measurements). We 
found that for ten of them the flux varied by less than 50 per cent, 
although more typically on the 10 — 15 per cent level. This de- 
gree of variability is not abnormal for NLSl on hourly time scales; 
therefore it seems that most members of the S-sample are in a typ- 
ical flux state (of course, this is based on limited observations). 
The XMM-Newton flux of IRAS 13349+2438 diminished by more 
than ~ 70 per cent since the two ASCA observations. During the 
BeppoSAX observation, IRAS 13349+2438 was another ~ 70 per 
cent dimmer then during the XMM-Newton pointing. Due to our 
strict criteria, IRAS 13349+2438 has not been included in our C- 
sample, but it does appear to possess some high-energy complexity 
(Longinotti et al. 2003). It could be that in even lower flux states, 
such as that observed with BeppoSAX, the spectrum would appear 
even more complicated. 

For three members of the C-sample, IH 0707^95, NGC 4051 
and PG 1211+143, it can be established from long-term, multi- 
mission light curves that the objects were in low-flux states dur- 
ing the XMM-Newton observations (see Leighly et al. (2002) for 
IH 0707^95, Uttley et al. (2003, 2004) for NGC 4051, and 
Yaqoob et al. (1994) for PG 1211+143^). Similarly, from a 10-day 
long ASCA observation (Dewangan et al. 2000) and comparison 
with the BeppoSAX observation (Comastri 2000) it can be shown 

^ The X-ray flux of PG 121 1+143 is also low compared to the Ginga ob- 
servations (Lawson & Turner 1997). 



that IRAS 13224—3809 was also in an X-ray low-flux state during 
the XMM-Newton observation. 

Data are sparse for PG 1535+547, PG 1402+261 and 
PHL 1092. Compared only to the ASCA fluxes, the XMM-Newton 
flux for PG 1535+547 was comparable while for PG 1402+261 it 
was ~ 40 per cent lower. For PHL 1092 the XMM-Newton flux was 
about 50 per cent lower than during the ASCA observation. It is dif- 
ficult to determine if this indicates a low-flux state for PHL 1092. 
Variations of this order are typical on hourly time scales in the 
0.1-2 keV range (Brandt et al. 1999; Gallo et al. 2004b); however 
the 2 — 10 keV variability is found to be negligible over similar 
time scales (Gallo et al. 2004b). 

If not convincing, it appears rather suggestive that the S- and 
C-samples defined in this work are, by coincidence, representative 
of NLS 1 in X-ray typical- and low-flux states, respectively. There- 
fore, it also follows that objects could transit from one sample to 
another depending on their X-ray flux state. 



4 OPTICAL AND UV PROPERTIES 

An advantage afforded with XMM-Newton is the availability of op- 
tical and UV data obtained with the Optical Monitor (CM; Ma- 
son et al. 2001) simultaneously with X-ray observations. For many 
of the NLSl observations, OM imaging data were also collected 
at UV or optical wavelengths. Therefore, whenever possible, we 
determined rest-frame 2500 A luminosities and optical-to-X-ray 
spectral indices (Uox)- 

The standard definition of Uox is: Uox ~ 
log{fx/ fuv)/log{vx/i^uv), where fx and /„„ are the intrin- 
sic flux densities at 2keV and 2500 A, respectively; and Ux and 
Uuv are the corresponding frequencies. In many of the observa- 
tions the rest-frame 2500 A was directly observed within one 
of the OM broad-band filters. Consequently, the flux density 
at 2500 A was estimated from the source count rate in that 
filter (Chen 2004). However, when this was not the case it 
was necessary to extrapolate from some measured UV wave- 
length (u) to 2500 A. The modified expression for aox is then: 
ctox = log{fx/ fu)/log{vx/vu) + log{fu/ fuv)/log{vu/vuv), 
where the second term on the right-hand side is the spectral 
slope between the measured UV flux density at u and 2500 A. 
If measured by Constantin & Shields (2003), the source-specific 
value for the UV spectral slope was used (the value shown in 
Table |2j. If the source-specific value was not known we adopted 
the spectral slope of the composite SDSS quasar spectrum between 
1300 - 5000 A {an = -0.44; Vanden Berk et al. 2001). 

The 2500 A flux density was measured from the available pre- 
processed pipeline products for each source. For most observations, 
multiple images were taken in the same filter, in which case the 
average source count rate was used to estimate the flux densities 
(Chen 2004). For Mrk 1044, Mrk 896 and Mrk 493, images were 
only made in the l/-band (5100 — 5800 A). For these observations 
an extrapolation to 2500 A was not done given that the optical spec- 
tra of AGN show a break in the continuum slope at about these 
wavelengths, which would introduce systematic uncertainties in the 
measurements of these three NLS 1 . Host-galaxy contribution was 
not excluded; thus will introduce a level of uncertainty in the re- 
ported values. 

In addition, FWHM(H/3) and Fe Il/H/3 ratios were collected 
from the literature for as many objects as possible. These are also 
included in Table|2| 
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Table 2. Optical and UV properties. The source name is given in column 1. The other parameters are the: 2500 A flux density (x 10 erg cm" 
(column 2); Fe II/H/3(column 3); FWHM of H/3 (kms"! ) (column 4); UV flux density (xlQ-^e ergcm-2 g-i Hz^^ ) (column 5); ~ 1100 - 4000 A 
spectral index (F oc u""^ ) (column 6); and ciox (column 7). Values in columns 5 and 6 are listed only when the 2500 A flux density (column 2) was not 
directly available from the data. 



General NLSl Sample (S-sample) 



(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


Source 


fuv 


Fe ll/H/3 


FWHM(H/3) /„ 


au 


aox 


Mrk 766 


0.57 


1.56 


1 100 - 




—0.958 


Mrk 359 


1.41 


0.50 


900 1.33 


—0.44 


— 1.276 


Mrk 1044 




0.77 


1310 — 






Akn 564 


1.31 


0.67 


865 - 




—0.941 


Mrk 896 




0.50 


1135 






Mrk 335 


5.43 


0.62 


1710 5.27 


-0.64 


-1.298 


NGC7158 






2100 






Mrk 493 




1.16 


800 






IZw 1 


1.61 


1.47 


1240 1.20 


-1.75 


-1.170 


Ton SI 80 


4.40 


0.90 


970 3.88 


-0.76 


-1.417 


PG 1448+273 


0.98 


0.94 


1330 




-1.407 


UGC 11763 


2.61 


0.63 


2210 




-1.489 


RXJ0323.2^931 


0.17 


0.65 


1680 




-1.131 


Mrk 478 


2.35 


0.97 


1630 




-1.495 


11 Zw 177 






1176 






IRAS 13349+2438 


1.06 


1.25 


2800 0.54 


-3.23 


-1.359 


PKS 0558-504 




1.60 


1250 






PG 1115+407 


0.82 


0.98 


1740 0.84 


-0.44 


-1.371 


Nab 0205+024 


1.02 


0.62 


1050 1.03 


-0.44 


-1.340 


PG 2233+134 


0.53 


0.89 


1740 




-1.407 


E 1346+266 


0.03 


0.98 


1840 




-1.048 


Complex NLSl Sample (C-sample) 


NGC 4051 


3.64 


0.94 


1170 




-1.484 


PG 1535+547 




0.47 


1480 






IH 0707^95 


1.54 


2.77 


1000 1.44 


-0.46 


-1.763 


IRAS 13224-3809 


0.65 


2.40 


650 0.60 


-0.44 


-1.536 


PG 1211+143 


2.89 


0.50 


1900 




-1.517 


PG 1402+261 


1.87 


1.10 


1623 




-1.441 


PHL 1092 


0.36 


1.81 


1300 0.30 


-0.44 


-1.560 



5 MULTIWAVELENGTH PARAMETERS 

In order to probe the multiwavelength behaviour of these NLSl 
we collected various parameters from the literature such as: radio 
loudness, optical emission line strengths, C IV strength, black hole 
mass and Eddington luminosity ratios. Comparing the two samples 
in all parameter spaces possible with these data reveled no clear 
difference between the C- and S-sample; thus these parameters will 
not be considered further. 



5.1 Contemporaneous X-ray and UV properties 

Measurements of aox from simultaneous UV and X-ray data 
clearly demonstrate a difference in the slope of the optical-X-ray 
continuum in the two samples (Figure|2j. The average index for the 
S-sample is -1.268 ± 0.179 compared to -1.550 ± 0.112 for the 
C-sample. A Kolmogorov-Smimov (KS) test comparing the two 
distributions yields a probability of < 0.1 per cent that they are 
drawn from the same sample. 

In Figure |2| we have plotted aox as a function of 2500 A 
monochromatic luminosity. For comparison, the UV luminos- 
ity dependence of aox for radio-quiet, type 1 AGN (aox = 



-0.136L„„ + 2.616) is also shown (Strateva et al. 2005). The S- 
sample appears to follow this relationship relatively well, but most 
of the C-sample fall below the Strateva relation. In combination 
with the fact that the C-sample objects are in X-ray low-fiux states 
(Section 3) this clearly demonstrates X-ray weakness in these com- 
plex objects during the observations. 

For five NLSl in the overall sample, at least two X-ray ob- 
servations with contemporaneous OM data are available. Three of 
these objects (Ton S 1 80, Mrk 478, Akn 564) are from the S-sample, 
the other two (IH 0707-495, NGC 4051) are from the C-sample. 
The epoch-specific aox and Luv were calculated for each of these 
observations. The results are plotted in Figure|2|as crosses. For the 
objects in the S-sample, the fluctuations in aox are small and per- 
haps even negligible if a complete error analysis was considered. 

In contrast, the aox variations are much more significant for 
the two objects in the C-sample (IH 0707^95 and NGC 4051). 
Figure |2| also indicates that the variations in aox are driven by 
changes in the X-ray flux, but that the UV flux also changes. For 
example, during the second observation IH 0707-495 was about 
six times brighter in the X-rays, while its UV flux diminished by 
~ 50 per cent. Increasing X-ray flux and decreasing UV flux was 
also seen in NGC 4051. 
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Figure 2. aox as a function of 2500 A luminosity measured from contem- 
poraneous X-ray and UV data. The filled (black) circles mark the C-sample 
NLSl and the open (red) squares are the S-sample NLSl. For comparison, 
the relation aox = —0.136Luu +2.616, found for radio-quiet type 1 AGN 
(Strateva et al. 2005) is shown as a green, solid line. The general NLSl sam- 
ple appears to follow the Strateva relation well, but the complex objects are 
X-ray weak. Multiple XMM-Newton observations exist for some of the ob- 
jects in our sample. These are shown as blue crosses. The Oox variations 
displayed by IH 0707-495 and NGC 4051 indicate significant changes in 
the X-ray flux and more moderate variations in the UV. Akn 564 appears 
excessively X-ray strong, but this is likely due to reddening of the UV spec- 
tram (Crenshaw et al. 2002). 

We also note that in the second observations of IH 0707^95 
and NGC 4051, both objects portray simpler high-energy spectra. 
That is, the null hypothesis for a power law plus Gaussian fit for 
IH 0707-495 and NGC 4051 at these epochs was 0.176 and 0.456, 
respectively. By our definition, neither object would have been con- 
sidered 'complex' at those times, suggesting that when measured in 
"typical" flux states, these source also follow the Strateva relation. 

5.2 Optical properties 

In Figure|3l the Fe Il/H/3 ratio is plotted against the FWHM(H/3). 
The average ratio for the S-sample is 0.93 ±0.34 with modest scat- 
ter as indicated by the reported standard deviation. The average of 
the C-sample is 1.43 ± 0.91 with considerable more scatter. The 
two samples clearly overlap. From a KS-test it is determined that 
the probability that they are drawn from the same sample is 8 per 
cent. Although not highly significant there are indications that the 
strongest Fe II emitters are the ones that also display high-energy 
spectral complexity. 



6 DISCUSSION 

6.1 X-ray weak NLSl 

According to Figure |2| objects that possess complexity in their 
2.5 — lOkeV spectrum appear to be X-ray weaker (i.e. they are 
found below the average Luv — Q-ox relation). At least in two of the 
objects (NGC 4051 and IH 0707-495) for which multiple, simul- 
taneous, X-ray/UV observations are available, extreme epoch-to- 
epoch changes are seen in aox , demonstrating that when the spec- 



Figure 3. The Fe ll/H/3 ratio is plotted against the FWHM(H/3). The sym- 
bols are as defined previously. The average ratio for the S-sample (open 
squares) and the C-sample (filled circles) are marked by a dashed and solid 
line, respectively. Although there is considerably more scatter in the aver- 
age ratio of the complex NLS 1 sample, the objects appear to possess much 
stronger Fe II emission compared to H/3. The average ratios are 0.93±0.34 
(general sample) and 1.43 ±0.91 (complex sample). The sample is limited, 
but it appears the most extreme Fe 11 emitters (e.g. Fe ll/H/3> 1.7) also dis- 
play complexity in their high-energy spectra. 

trum appears less complex. X-ray emission is stronger and olox is 
more typical. 

The transiting of objects from the C- to the S-sample with in- 
creasing flux is predicted by both partial-covering and reflection 
models. In terms of partial-covering the increased spectral com- 
plexity during the low-flux states occurs because the intrinsic spec- 
trum (presumably a blackbody plus power law) is highly absorbed, 
imposing edges and curvature in the detected spectrum. For reflec- 
tion, the low-flux spectrum will likely be reflection dominated and 
the spectral curvature and drops will be associated with various re- 
flection components (predominately the blurred iron line and soft 
excess; e.g. Ross & Fabian 2005). 

6.1.1 Absorption 

The X-ray weakness is, of course, in agreement with absorption 
models such as partial-covering where the X-ray emitting region is 
more absorbed than the UV region. This could occur, for example, 
if the absorber is located somewhere between the X-ray and UV 
emitting regions of the accretion disc. Consequently, UV fluctua- 
tions could arise from e.g. changes in the accretion rate, while the 
X-rays variations will appear larger because of the addition affect 
of the changing absorption. Therefore the extreme changes in the 
X-ray flux due to absorption can drive the fluctuations in aox when 
compared to the relatively unabsorbed and less variable UV flux. 
We also note that in the multi-epoch observations of IH 0707-^95 
and NGC 4051 the UV flux diminished when the X-ray flux was 
higher. 

A prediction of the partial-covering model (e.g. Tanaka et al. 
2004) is that changes in the covering fraction or column density 
of the absorber(s) will account for, at least, the long-term spectral 
variations; but that the unabsorbed, intrinsic flux will remain rel- 
atively constant over time. We can test this prediction by reexam- 
ining the two XMM-Newton observations of IH 0707^95 (Gallo 
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et al. 2004a). During both observations, the absorbed component 
of the intrinsic spectrum was seen through a column density of 
^ 10^^ cm~^ ; however during the first, low-flux state observa- 
tion the covering fraction^ of the absorber was much greater. The 
intrinsic, 2 keV flux density during the first and second observation 
of IH 0707^95 was 13.13 and 8.87 x 10"^° ergcm'^ s'^ Hz'^ , 
respectively. This yields aox ~ —1.17 at both epochs, demonstrat- 
ing that a common intrinsic flux modified by partial obscuration is 
possible in accounting for the differences in brightness over time. 



6.1.2 Reflection 

Steep and fluctuating aox is not inconsistent with models, such as 
light bending (e.g. Miniutti & Fabian 2004), which incorporate gen- 
eral relativistic effects close to the black hole to account for some of 
the X-ray variability. In these models the low-flux, reflection domi- 
nated spectrum results when most of the emission from the primary, 
power law component is bent back toward the black hole and never 
reaches the observer. In comparison with the relatively constant UV 
emission (most likely produced by the colder accretion disc) aox 
will be smaller (X-ray weaker). As the emission from the power law 
continuum becomes important, the X-ray source becomes brighter 
and aox begins to flatten. 

6.2 Optical Fe 1 1 emission 

Not all of the NLSl in the C-sample possess strong optical Fe II 
emission (specifically Fe Il/H/3), but it does appear that the most 
extreme Fe II emitters (e.g. Fe II/H/?^ 1.7; IRAS 13224-3809, 
IH 0707-495, PHL 1092) are spectrally complex. 

Firm conclusions cannot be drawn due to the small number 
statistics, but this raises the possible conjecture that not all object 
which display complex X-ray spectra possess strong iron emission, 
but possibly all objects with extreme Fe II emission do possess 
the mechanism to create high-energy complexities. A simple test 
would be to make X-ray observations of a few extreme Fe II emit- 
ters, say with Fe Il/H/3> 1.7, as suggested by Figure|3| and deter- 
mine if their high-energy spectra are complex. 

However, the fact that a source can transit from one sample to 
another depending on its X-ray flux state (Section 3) raises further 
issues and limits the conclusions we can make regarding Fe II emis- 
sion. The primary issue is whether the Fe II emission is variable 
and if it depends on the X-ray continuum flux. There are several 
studies which show evidence for (e.g. Kollatschny & Fricke 1985; 
Kollatschny et al. 2000; Vestergaard & Peterson 2005; Wang, Wei 
& He 2005) and against (Goad et al. 1999; Kollatschny & Welsh 
2001) significant Fe II variability. Since the Fe II measurements 
presented here are not contemporaneous with the X-ray measure- 
ments, we can only speculate on a possible connection between 
Fe II strength and X-ray spectral complexity. 

An alternative possibility is that the strongest Fe II emitters 
simply show the greatest variability. If so, this makes it more likely 
to catch the strongest Fe II emitters in an extreme (low) flux state, 
as opposed to an average NLSl. IRAS 13224-3809, IH 0707^95 
and PHL 1092 are certainly know for there strong Fe II emission 
and significant variability. From the current analysis it is difficult to 
distinguish the two theory. 

^ Note that the covering fraction and column density are degenerate in 
partial-covering models and cannot be distinguished. 



7 SUMMARY 

We have investigated the nature of NLS 1 which exhibit complexity 
in their high-energy (2.5 — 10 keV) spectrum. We identified seven 
out of twenty-eight NLSl whose high-energy spectrum show sig- 
nificant deviations from a simple power law. We compared multi- 
wavelength properties of this complex sample of seven NLSl with 
the larger sample to investigate whether we could identify any un- 
derlying physical processes that could be responsible for the spec- 
tral complexity. The main results follow. 

• Considering long-term X-ray variability of the sources, the 
complex sample of NLS 1 (C-sample) seems representative of ob- 
jects in X-ray low-flux states, whereas the NLSl in the general sam- 
ple (S-sample) appear to be in a typical flux state. 

• In cases where multi-epoch measurements of aox were pos- 
sible, it appeared that aox approached more typical values as the 
complexity in the high-energy spectrum diminished. For IH 0707- 
495 and NGC 405 1 (both in the complex sample) the variations in 
aox were extreme and while predominately due to variability in the 
X-rays, the UV flux did also change. Specifically, in these two ob- 
jects, the UV flux actually diminished when the X-ray flux was high 
and the spectral complexity less prominent. This is not inconsistent 
with either partial-covering or reflection scenarios. 

• On average, the C-sample showed stronger Fe II emission 
(Fe Il/H/3) than the general sample. Specifically, the three most 
extreme Fe II emitters (with Fe Il/H/3> 1.8) of the twenty-eight 
NLSl, were all part of the C-sample of only seven. This raises an 
intriguing possibility that we may be able to identify complicated 
X-ray NLSl based on the strength of the optical Fe II emission. 
It is not certain from this current analysis if the possible connec- 
tion between X-ray spectral complexity and Fe II strength is due 
to high iron abundances or because strong Fe II emitters may ex- 
hibit greater X-ray variations and consequently are more likely to 
be caught in an extreme (low) flux state. 

Both absorption and reflection models predict the observed 
behaviour in aox- It is interesting to note that in the case of 
IH 0707^95 a common, intrinsic aox is predicted by the partial- 
covering model and is consistent with the XMM-Newton observa- 
tions. 

Partial-covering models make no strong assumption on the na- 
ture of the intrinsic spectrum, only that it is partly obscured. In 
this sense reflection models, which can describe the broad-band 
X-ray properties are attractive. While it seems obvious that reflec- 
tion close to the black hole should occur it also seems premature 
to abandon concepts which employ absorption in the black hole 
environment to explain some of the X-ray behaviour. 
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